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Very often, when studying topological or dynamical properties of random scale-free networks, it 
is tacitly assumed that degree-degree correlations are not present. However, simple constraints, 
such as the absence of multiple edges and self-loops, can give rise to intrinsic correlations in these 
structures. In the same way that Fermionic correlations in thermodynamic systems are relevant 
only in the limit of low temperature, the intrinsic correlations in scale-free networks are relevant 
only when the extreme values for the degrees grow faster than the square-root of the network size. 
In this situation, these correlations can significantly affect the dependence of the average degree of 
the nearest neighbors of a given vertice on this vertices’s degree. Here, we introduce an analytical 
approach that is capable to predict the functional form of this property. Moreover, our results 
indicate that random scale-free networks models are not self-averaging, that is, the second moment 
of their degree distribution may vary orders of magnitude among different realizations. Finally, 
we argue that the intrinsic correlations investigated here may have profound impact on the critical 
properties of random scale-free networks. 
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I. INTRODUCTION 

Very often, the network representation of a real system 
displays the so called scale-free property M. These net¬ 
works typically exhibit a degree distribution in the form 
of a power law, P{k) ~ k~'^, with the most connected 
nodes reaching a degree k orders of magnitude higher 
than the average degree (fc). In fact, in most models 
of scale-free networks the degree of the most connected 
nodes diverges at the thermodynamic limit, K^ax ^ 
with N being the number of nodes in the network and 
0 some positive exponent. If 7 < 3 and with K^ax di¬ 
verging, the second momenh (/c^), also diverges, leading 
to some surprising results [^, such as the critical point 
of several models in scale-free networks going to zero in 
the thermodynamic limit im. 

The exponent 9 controlling the divergence of Kmax can 
also affect signihcantly the structure of the network [T^ . 
Simple extreme value statistics [l3l - [T5j shows that, even 
if no other constraint is imposed, the most connected 
node in the network should have a degree of the order 
of Kjnax ~ iV®", with 9n = 1/(7 — 1 ), which represents 
the natural cutoff of scale-free networks [l^. However, 
using the natural cutoff may lead to unexpected prop¬ 
erties in the network structure. Typically, if the links 
are placed at random, one should expect the average 
number of links between a pair of nodes do be given by 
nij = kikj/N{k) [il[i3- If the nodes i and j are ex¬ 
tremely connected nodes, one may expect their degree 
to be of the order of K^ax- In this case, for scale-free 
networks with 7 < 3, the maximum value for grows 
with the network size and should reach values above one. 
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that is, for these scale-free networks one should expect 
multiple connections between the most connected nodes. 
However, it is usually the case that network models do 
not account for multiple connections, which would lead 
to the presence of intrinsic correlations [l^ . As a matter 
of fact, these intrinsic correlations have been observed 
in different random network models [Tsl [T^. Moreover, 
since they can be related to correlations of Fermionic sys¬ 
tems due to the exclusion princi ple, a practical analyt¬ 
ical treatment becomes possible [Ol- More specifically, 
in analogy with correlations of Fermionic systems that 
may become irrelevant at sufficiently high temperatures, 
the intrinsic correlations in networks may practically dis¬ 
appear if a more restrictive constraint is imposed to the 
extremely connected nodes, Kmax ^ with 9s < 1 / 2 . 
In this case, the maximum expected value for riij does 
not diverge with N and multiple connections would not 
be present even if allowed. This so-called structural cut¬ 
off [il] represents a constraint that is strong enough to 
turn degree-degree correlations in the network negligible. 

In this work we investigate the problem of intrinsic 
correlations in random scale-free networks. This is per¬ 
formed by measuring the dependence of the average con¬ 
nectivity of the nearest neighbors with the degree of the 
node, knn{k)- We confirm that for scale-free networks 
the intrinsic correlations lead to a disassortative network, 
with knn being depreciated in the most connected nodes. 
We observe, however, that the particular form for 
varies widely depending on the network instance, that 
is, two networks of same size, generated with the same 
underlying degree distribution, may display highly dif¬ 
ferent dependencies for fc„„(fc). From our results, we ob¬ 
serve that the moments of the degree distribution, (k) 
and (fc^), are determinant in the shape of knn{k). We 
then use an analytic approach, based on the method of 
Fermionic correlations to obtain the functional de- 
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pendence of knnik) given the network size and degree 
distribution. Since the moments of the degree distribu¬ 
tion vary from one realization of the network to another, 
we need to adjust the minimum and maximum cutoffs 
of the degree distribution to each realization. Using the 
properly adjusted degree distribution we observe that our 
approach is capable of correctly predicting the functional 
form of knn{k). Our results show that random scale-free 
network models are not self-averaging. More precisely, 
relevant parameters of the network, such as the second 
moment of the degree distribution, can vary orders of 
magnitude for different realizations. 


II. THE INTRINSIC CORRELATIONS 


To build our scale-free networks we use the configura¬ 
tion model (l9l - [^ . In this model, we start by choosing 
the degrees of each of node from a power-law distribu¬ 
tion, namely, we let the degrees of the nodes to be given 
by, 


ki = 


[( 


^- 7+1 _ r 

max min J ^ 


K. 


- 7-1 


( 1 ) 


where [-J is the floor function and are random numbers 
uniformly distributed between 0 < < I. One should 

note that, while the degrees ki are integers, the parame¬ 
ters Kmax and Kmin may assume arbitrary values. After 
the degrees are determined, the network is constructed 
by randomly selecting a pair of nodes. The probability of 
selecting a particular node should be proportional to its 
remaining degree. Once the pair is selected, a connection 
should be placed between them and, if the connection is 
allowed, their remaining degree is decreased by one. Note 
that we include two important constraints, namely, there 
is no connection between a node and itself, and there 
is no redundant connections, that is, each pair of nodes 
will connect at most once. This later constraint is simi¬ 
lar to the exclusion principle of Fermionic particles, being 
particularly important in scale-free networks with 7 < 3 
as it leads to the so-called intrinsic correlations. The 
imposed constraints may result in frustrated attempts 
to construct the network because at some point, while 
constructing the network, one may reach a situation in 
which all the nodes with remaining degree are already 
connected. In this case, this attempt is rejected and the 
whole process needs to be restarted by drawing a new set 
of degrees. 

After the network is built, we proceed to calculate the 
average nearest neighbor degree This is done by 

computing the average degree of the nearest neighbors of 
each nodes, and then averaging over nodes that possess 
the same degree. In Fig. [1] we show the results obtained 
for two network realizations generated with the same set 
of parameters, number of nodes N = 10®, 7 = 2.5, and 
Kmin = 2. We tune the parameter Kmax+ = N — 1^ how¬ 
ever, since frustrated attempts are discarded, the most 
connected nodes of the network never reach this limit. 
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FIG. 1: Average connectivity of the nearest neighbors of ver¬ 
tices with degree k, considering two values of (k^) and the 
same value for the mean degree (k). The symbols represent 
simulations considering the configuration model . The 

dashed lines represent the results of the analytic approach de¬ 
scribed in the text. Here we considered N — 10® for the size 
of the network and 7 = 2.5 for the degree distribution. 


Curiously, these two realizations of the model results in 
quite distinct forms for although both present 

a clear disassortative trend, with knn decreasing for the 
most connected nodes. The difference between these two 
realizations can be elucidated by noting that, although 
having similar average degrees, they present very differ¬ 
ent values for the second moment of the distribution, 

(k^)- 

In Fig. [ 2 ] we show another two realizations obtained 
with the same conditions, but now we selected realiza¬ 
tions with similar second moments but distinct av¬ 
erages (k). The results for again do not match in 
this case, although the difference is not as pronounced as 
in Fig. [TJ The fact that the results vary from realization 
to realization indicate that the network model is not self¬ 
averaging. In this way, to understand the effect of the 
intrinsic correlation, we need more information than just 
the network size and the value of the exponent 7 . The 
absence of the self-averaging property can be observed in 
Fig. ([3]), where the probability distributions for the mo¬ 
ments (k) and (fc^) are shown for a set of 10 ® samples of 
networks built from the configuration model and consid¬ 
ering the Fermionic conditions. It is important to notice 
that the range of the (fc^) distribution covers two orders 
of magnitude. 


III. THE SELF-CONSISTENT EQUATIONS 

By considering that a network is uncorrelated results 
in an average number of connections between a pair of 
nodes i and j given by = kikj/N{k). As mentioned 








3 



10 ^ 

k 


FIG. 2: Average connectivity of the nearest neighbors of ver¬ 
tices with degree k, considering the same value for cumu- 
lant and two different values for mean degree {k). The 

symbols represent simulations considering the configuration 
model [l9l-[^. The dashed lines represent the results of the 
analytic approach described in the text. Here we considered 
N = 10® for the size of the network and 7 = 2.5 for the degree 
distribution. 
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FIG. 4: The g{k) functions obtained from the self-consistent 
set of equations 3iS C3,lcill&ti0ci for I^min — 2 BjIIcI I^max — 
Here the following cases have been considered: 7 = 
2.5 (circles), 7 = 2.6 (rectangles), 7 = 2.7 (triangles), and 
7 = 2.8 (stars). 

can be written as. 



FIG. 3: The probability distributions for the moments {k) 
and (fc^). Here we generate 10® samples of networks from the 
configuration model, but do not allow for self and multiple 
connections. 


before, this relation cannot hold for scale-free networks 
with 7 < 3, unless some restriction is imposed to the 
maximum degree of the nodes, Kmax- Our goal here 
is to be able to determine, given a network with distri¬ 
bution P[k) nodes of degree fc, the dependence on k of 
the average nearest neighbors The effect of intrin¬ 
sic correlations on random scale-free networks has been 
previously studied in [13 , where it was shown that the 
average number of connections between a pair of nodes 


^ 9{ki)9{kj) 

^ + 9{ki)9{kj)' 


( 2 ) 


where the form of the function g{k) has to be determined 
by the consistency of ki = From that, one ob¬ 

tains the following self-consistent set of equations for the 
values of g{k): 


g{k) = k 



P{k)g{k') 
g{k')g{k) -b 1 


-1 


( 3 ) 


In order to obtain g{k), we translate Eq. ([3]) in terms of 
an iterative process. We start with go{k) = k/{N 
as an tentative guess for g{k)^ and place go{k) on the right 
side of Eq. ([3]) to obtain gi[k) on the left side. We then 
proceed to obtain g 2 {k) from gi{k) and so on, until we a 
self-consistent solution is obtained. In Eig.|3]we show the 
form of g{k) obtained through this process for different 
expressions of the degree distribution P{k). Erom g{k), 
it is possible to obtain the dependence of the average 
nearest-neighbor degree as, 


knn{k) 


9{k) 

k 




g{k')k’ 

g{k'){g{k) + !)■ 


( 4 ) 


However, as we observed from the results for the configu¬ 
ration model, two important factors need to be taken into 
account. First, some attempts to build the networks are 
rejected, therefore the effective degree distribution P{k) 
may be distorted by this selection criteria. Second, since 
the model is not self-averaging, the results for a given 
network may differ greatly from the average result given 












4 


the degree distribution used to generate the network. To 
obtain we need to determine the dependence we 

will introduce in the degree distribution P{k) to match 
each network realization. The simplest approach is to 
set P{k) exactly as the number of nodes with degree k 
that we obtained in a given realization. Our simulations 
show that the obtained form for knn{k) follows closely 
the results observed with the configuration model (data 
not shown). 

Alternatively, instead of matching each value of P{k) 
we tried to reduce the number of parameters by matching 
only the first and second moments of the degree distribu¬ 
tion. To do that we observe that from Eq. m we obtain 
that, on average, 


- (k + T)“'>'+^ 


k: 


J-^max 


( 5 ) 


However, considering that Kmin can be non-integer, the 
minimum degree can be written as kmin = [Kmin\ ^ and 

p{km.n) = N — — ■ (s) 

P-min P-max 


Therefore, by varying continuously the parameter Kmin, 
we reduce continuously the expected number of nodes 
with the minimum degree P{kmin)- Now, given the num¬ 
ber of nodes N and the exponent 7 , the two free param¬ 
eters, Kmin and Kmax, Can be set to match the values of 
the first and second moments obtained from a particular 
random network realization. The results from the con¬ 
figuration model presented in Figs. © and @ are com¬ 
pared with the expected forms for the average nearest- 
neighbor degree obtained with this approach. As 


depicted, there is good agreement between the distinct 
approaches. 


IV. CONCLUSION 


In this paper we have studied the intrinsic correlations 
of scale-free networks considering a self-consistent set of 
equations g{k), in order to analyse the dependence of the 
average connectivity of nearest neighbors with the degree 
of nodes. We have confirmed the disassortative behavior 
of scale-free networks, and demonstrated that the par¬ 
ticular form for is strongly dependent on the two 

first moments of the degree distribution. Moreover, us¬ 
ing the properly adjusted degree distributions, based on 
the minimum and maximum cutoffs of the original degree 
distribution, we have observed that out approach is capa¬ 
ble to provide correct predictions for the functional form 
of knn{k). Finally, our results also indicate that random 
scale-free models are not self-averaging since these mod¬ 
els are highly determinant parameters of the network, 
such as the second moment of the degree distribution, 
may vary orders of magnitude from one realization to 
the next. Therefore, it is likely that our results have sig¬ 
nificant implications on the critical behavior of models 
defined on complex networks substrates [H, [23 - [^ . 
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